Nonlinear resonances in 5-kicked Bose-Einstein condensates 
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We investigate the Quantum Resonance (QR) regime of a periodically kicked atomic Bose-Einstein 
condensate. We find that the clearest indicator of the nonlinear dynamics is a surprisingly abrupt 
cut-off which appears on the main "Talbot time" QR. We show that this is due and excitation 
path combining both Beliaev and Landau processes, with some analogies to nonlinear self-trapping. 
Investigation of dynamical instability reveals further symptoms of nonlinearity such as a regime of 
exponential oscillations. 
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The interaction between nonlinear dynamics and quan- 
tum dynamics has stimulated new domains of cold atom 
physics. For instance, the suppression of chaotic classical 
diffusion in a corresponding quantum system evolving un- 
der the linear Schrodinger equation has been investigated 
experimentally using cold atoms subjected to regularly 
spaced short pulses ((5-kicks) from standing waves. But 
studies of the quantum (5-kicked systems, which have as a 
classical analogue the famous nonlinear paradigm which 
is the Standard Map, had certain limitations: probing 
the dynamics close to the chaotic classical limit requires 
small values of Planck's constant, h. Even though h in 
the cold atom experiments is a scaled (and hence vari- 
able) effective value, the h ^ regimes remained out of 
reach. However, an important advance was the realiza- 
tion in 2002 [H that atoms kicked with periods which are 
a rational multiple of the Talbot time T = An (so-named 
because of an analogy to the Talbot optical effect) repre- 
sent a somewhat pathological case. They correspond to 
a classical 'image' system where the effective value of h is 
played by the dephasing from resonance and so can take 
very small values. This has stimulated a series of deli- 
cate experiments by several groups worldwide probing 
the rich physics of both the Quantum Resonance (QR) 
regime as well as the closely related Quantum Acceler- 
ator (QA) regime (which is in effect, quantum resonant 
behavior in the presence of gravitational acceleration). 
Further theory [3] includes proposed applications such as 
the realization of a quantum random walk algorithm 

Since well-defined initial momenta are needed, most 
recent QR or QA experiments employed i5-kicked Bose- 
Einstein condensates (BECs). But this in turn opens up 
a new and quite different possibility: namely the largely 
unexplored regime where nonlinear dynamics, arising 
from the many-body nature of the BEC, combine with 
the (5-kicked quantum dynamics. To date, the deep un- 
derstanding of collective excitations, acquired from other 
areas of BEC physics, has not been applied to the unique 
dynamical features of the (5-kicked atoms. 

A few theoretical studies have considered interactions 



in QR/QA: in ^5j] dephasing of resonant transport at the 
Talbot time as a function of nonlinearity parameter g 
was found; in Q exponential growth of non-condensate 
atoms was predicted for certain parameters at half the 
Talbot time; in [3] significant differences were reported 
for QA dynamics between attractive and repulsive inter- 
actions. Only recently, though, was it demonstrated Q 
that an approach based on Bogoliubov phonon modes 
is essential; in particular, the onset of dynamical insta- 
bility was attributed to parametric resonances. Para- 
metric instabilities through periodic driving of collective 
modes have been investigated in several BEC studies f9| . 
Higher-order effects, resulting from phonon-phonon in- 
teractions, have been experimentally studied, e.g., in the 
context of excitation lifetimes and in a nonlinear cou- 
pling between two phonon modes [l^ . Beliaev and Lan- 
dau (BL) couplings provide the dominant contribution to 
such nonlinear mode conversions; their importance was 
demonstrated in a series of recent experiments with BECs 
excited by optical Bragg pulses [3] ■ 

Here we have, for the first time, quantitatively investi- 
gated an interacting (5-kicked BEC by mapping the posi- 
tion and properties of the resonances of all its low-lying 
modes. The leading Talbot time resonances, which ap- 
pear to dephase and disappear at T = 47r, in fact are 
shifted and broadened (they are strong over AT ~ 1). 
But our key finding is that they acquire an extraordi- 
narily sharp cut-off at their maximum (AT 10"'^). We 
show this behavior is due to a nonlinear feedback process, 
originating from the resonant kicking and a combination 
of both Beliaev and Landau coupling. In other nonlinear 
resonances, we find novel features, not previously seen in 
other driven BECs, such as exponential oscillations and 
Fano-like profiles. We calculate the local Lyapunov ex- 
ponents and model the Fano profiles quantitatively. The 
kicked-BEC experiments to date, have used effective val- 
ues of the nonlinearity (g ^ 0.5) which are only slightly 
smaller than those {g > 1) needed for the effects we find. 

For non-interacting single-particle dynamics, an at- 
tractive features of the (5-kicked system, for quantum 
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chaologists, is the extensively-investigated simple quan- 
tum map which stroboscobically evolves the system from 
kick n to kick n+l. For example, expressing our quantum 
state in a momentum basis, tpi^i t) — X]/^o 1^ 
write: 



a((n + l)r) = Ug=o(T) a(t = nT) , 



(1) 



where a(i = nT) is a vector with the amplitudes ai. 
The atoms experience a kicking potential Vkickix^t) — 
K cosxJ2n^(^ ~ nT). The corresponding unitary time 
evolution operator factors (exactly) into a free-evolution 
part Ufree and a kick part Ukick, i-e.. 



Ug=o(T) = Ufree Ukick 



(2) 



The first exponential term represents the free evolution 
under the kinetic energy operator in some units where 
the atom mass M = 1 and h = 1. Clearly, since the 
atomic momentum / — 0,±1,±2... is quantized in units 
of the recoil momentum, T = 47r implies U free = 1 for all 
momentum states, so consecutive kicks add in phase. The 
result is a phase-matched absorption of energy from the 
field yielding ballistic transport. In contrast, the non- 
resonantly kicked atoms experience diffusive growth in 
energy. 

We consider a uniform, tightly-confined, effectively ID 
BEC with periodic boundary conditions. The corre- 
sponding field-free, many-body Hamiltonian, in a mo- 
mentum representation, is: 
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9iD sr^ .f-t- - fo\ 

2^ "■l%^r,ia:j+l-ra , (3) 



where = h^k'^/{2M) and L is the BEC size. The 
ID interaction constant gm — 2hasUj±, depends on 
the atomic scattering length as and the transverse trap 
frequency luj_. In the presence of linearized perturba- 
tions around the macroscopically occupied k = ground 
state, H may be diagonalized up to quadratic order by 
the Bogoliubov transformation dk ~ u^bk — Vf;P_i,, for 
k ^ 0, with Uk — Vk = l/{uk + Vk) = y/huJk, where 
huJk — [ek{ek + 2giDn)]^^^ and n — N/L is the atom den- 
sity. Then we can expand H in the orders of ^/N : H = 
const + £r(2) _^ ^(3) _^ ^(4) ^T^^^^ fj{2) ^ Y^^_^^ hwkblbk 

and the cubic part H^^^ describes the leading order con- 
tribution to the interactions between phonons 

^ ^ ^(r,p6t St 5t_^_^ + /^Jlb\k+, + H.c.) , (4) 
q,p 

where q,p,{q + p) ^ and k, = y^guj/ L. The coef- 
ficient Tqp represents a process in which three phonons 
are created or annihilated and Aqp a process in which one 
phonon with momentum q + p decays into two phonons 
with momenta q and p (Beliaev term) or its inverse in 
which two merge to produce a third (Landau term) . Since 



the energy of the excitations around the ground state is 
positive, the processes described by Tgp are suppressed 
by energy conservation. In terms of the Bogoliubov am- 
plitudes we obtain 

Aqp — UpUqliq^p -)- 'ZUpVqVqJ^p 'ZUqVpUq^p VpVqVqJ^p 
Pgp — UqVpVqJ^p UpUqVqJ^p (5) 

We assume macroscopic occupancy for the low-lying 
modes of interest and that phase fluctuations of the ID 
condensate may be neglected. Hence we will treat the 
Bogoliubov mode amplitudes classically. The (5-kicked 
map ([2]), in the presence of the quadratic Hamiltonian 
i?^^^ alone, requires only a straightforward modification 
to its free-evolution part: 



\Jg{T)=B-^e-"^^ IBUk 



ick 



(6) 



Where B denotes the Bogoliubov transform, exp (— iojT) 
is a row vector with the mode frequencies and X is 
the identity. The eigenvalues of the non-unitary matrix 
Ug(r) indicate dynamical stability. It is easy to prove 
they in general come in quartets A, I/A, A*, I/A*. Then 
|Amaa:| > I (whcrc Xrnax is the largest eigenvalue and 
the local Lyapunov exponent) imply dynamical instabil- 
ity and exponential growth in the relevant modes (at least 
for short times). 

We compare Eq. ^ to the numerical solutions of the 
ID Gross-Pitaevski equation (GPE) (here rescaled to di- 
mensionless units fisj) for an initially uniform BEC: 



idt'4){x, t) = 



¥1 



■ g\%lj{x,t)\^ + Kcoax Ft] i^{x,t) 

(7) 

where Ft — 5{t — nT). We stu dy q ~ — 10, realistic 
with the current experiments 1^ 16l |. for X < I, which 
allows only the creation of discrete low-lying phonon ex- 
citations. 

Figllja) shows GPE numerics; it maps the average 
BEC response, < I - |ao(t)p >*= \ ELi 1 " \ao{nT)\^ 
averaged over the first t kicks, for K — 0.5. At g ~ 0, 
the Talbot time resonance at T = 47r is perfectly sym- 
metrical, (as are the fractional resonances on either side). 
For g ~ — I an asymmetry develops, due here to the 
lifting of the degeneracy between the lowest modes: the 
main resonance splits into a mode I resonance wiT ~ 2-k 
and mode 2 resonance IJJ2T ~ Stt. The mode 2 resonance 
rapidly decays away as the gap uj2 — oJi increases: direct 
couphng < 0|?7fcicfc|2 J^iK) with the condensate is 
small. A slight asymmetry was noted in the GPE numer- 
ics in Ref. [,5| for the Talbot time resonance at g ~ O.I, 
which we now attribute to this regime. 

But the most striking feature is the very sharp 'cut-off 
appearing at gr > I (and still exists even at 17 = 20). It is 
also evident in the second harmonic of the resonance (up- 
per half of the graph) . Fig[Hb) shows that even a grid of 
100 points per unit of T (each circle representing a GPE 
simulation for 30 kicks and g = 5) is too coarse to resolve 
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FIG. 1: (a) Losses from the unperturbed condensate < 1 — 
|ao(i)P >t averaged over time, calculated from the GPE, for 
kick strength K — 0.5. The BEC resonance evolves from the 
well-known Talbot time resonance of the QKR (at T — AttM, 
M — 1,2. . ., for g — 0) and develops an abrupt cut-off for 
any g > 1. (b)Illustrates the sharpness of the cut-off. GPE 
calculations with even 100 points (black stars) for every unit 
of T still fails to resolve the cut-off. Plot for g = 5, K = 0.1. 
The quantum map excluding BL processes, incorrectly gives 
a symmetric, unshifted resonance. Inclusion of BL coupling 
perfectly reproduces not only the shift but also the cut-off; the 
simple Eq. ([UJ also provides reasonable agreement (inset). 



the cut-off (in comparison, such a grid could resolve the 
famously narrow g = Talbot-time resonance). 

The dotted line shows the Bogoliubov map ^ here 
incorrectly produces a symmetric resonance and fails to 
shift the resonance away from loiT = 27r. Hence we 
need to include the neglected phonon-phonon interac- 
tion terms between the kicks from H^^K The free-ringing 
exp (— iwT) part of the map [Eq. ^] must be replaced 
by a set coupled equations following from the Heisen- 
berg equations dtbk — —iujkbk — i[bk, H'^^^ /H], for k ^ 0, 
where we replace bk by the rescaled classical amplitudes 
bk {bk)/VN. Figure [TJb) shows only four lowest ex- 
citations (fc = ±2, ±1) provide an excellent quantitative 
agreement with the GPE, accurately reproducing the cut- 
off for small depletion (< 10%) of the ground state. 

We simplify further by transforming (by symmetry 
bk = &-fe) to the basis |/) + \ - I)), for / 7^ 0, so 

that {l\Uk^ck\n) = Uni = z'-"J„_/(if) +z-('+")j,+„(X) 
if n,l> 0, but Uoi = V2i-^Ji{K) and Uqo = MK). We 
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FIG. 2: Upper figure (left): Probability for mode 2 averaged 
over 100 kicks for K = 0.5; (right) Shows |Ama2,|, largest 
eigenvalue of Ug. Bright regions denote |A| > 1 and hence 
exponential behavior (dynamical instability). The unstable 
Ni, N2, ..Nn series of nonlinear resonances (which only ap- 
pear for g > 1) correspond to {ui + UJ2) — 2miT. The as- 
terisk denotes position of the 'instability border' found by 
which is thus due to The L series are resonances 

which evolved from partial or full resonances of the Talbot 
time at g = 0. They are stable, but in spite of this are much 
stronger than the exponential resonances. Lower figure (left). 
Mode 2 probability of N^, near the asterisk of upper figure. 
T = 6.12, g = 2.5. The exponential growth persists for only 
a finite time; it is then replaced by exponential decay, leading 
to exponential oscillations{\og scale shown in inset) . Lower 
figure(right). Mode 2 near The cluster of 3 interacting 
resonances have "Fano-like" profiles. These are well repro- 
duced by Eq. corrected with J/'-^^ in Eq. (g)), including 
only the 7 lowest modes. 



then only need two coupled equations to accurately re- 
produce the cut-off. Moreover, neglecting the small Tgp 
terms, we obtain (from the free-ringing plus BL terms) 

dth ^ -t[uji + 2CiRe(b2)]bi + 2C2bfb2, 

dtb2 = -tL02b2~t[Ci\h\^ + C2bl], (8) 

where Ci = k(A_i^2 + A2,-i), C2 = SAi 1, and R = 
g/{2y/2TT) (note A_i,2 = Ai _2, etc.). We can even 
obtain a reasonable cut-off if we set the smaller term 
C2 = 0. Hence, the main effect of mode 2 is simply 
to provide a phase-shift on uji . If we integrate &2 while 
keeping constant, we can replace the full map by: 

dtbi ~ -iujibi + iAsin^{LU2t/2)\bi\'^bi - ^7^(^) , (9) 

where A = AC^/lu2, n{t) = V2Ji{K){ui - vi)F{t). 
Figdjb) inset shows that we still get reasonable agree- 
ment with the full model for weak K = OA. In this 
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regime |6o(OP ^ ^ - ^(OP- Writing 61 = pe*^, a 
phase space analysis in the p, plane reveals a separa- 
trix curve which appears at the cut-off parameters and 
bounds the value of p. We may also describe the physical 
mechanism thus: if the — *■ 1 transition is initially only 
slightly off-resonant, the kicking starts populating mode 
1 effectively; with the increasing mode 1 population, BL 
processes begin to populate the empty modes / — ±2; 
this mode 2 population provides a phase-shift bringing 
mode 1 closer into resonance; the nonlinear feedback ac- 
celerates the growth in bi which in turns brings mode 

1 further into resonance. However, if the 0^1 tran- 
sition is initially too far off-resonant (beyond the sharp 
cut-off value) the nonlinear feedback cycle cannot start 
and the population oscillations between the modes are 
suppressed. An analogous model of two-mode dynamics 
with continuous driving rather than S kicks is reminis- 
cent of a macroscopic self-trapping effect in a BEC in a 
double- well potential [l3|, but in the present work the 
cut-off is considerably sharper. 

Figl2] maps the average probability of mode 2 (av- 
eraged over 100 kicks) for K = 0.5. The right hand 
side maps regions of dynamical instability |Amaa;| > 1- 
We analyze dynamical stability by mapping the eigen- 
values of Ug(T) for all the resonances of the lowest 3 
excited modes. We divide the resonances into (1) the 
"linear" family L{n, I) (ie those which evolve from the 
linear case and converge at g = to a rational frac- 
tion of the Talbot time. The resonance in Figllja) is 
the L{1,1) (first resonance of mode I — 1). (2) The 
"nonlinear" resonances Nn and Vn which vanish in the 
absence of interactions, at g = 0; the Nn correspond to 
(wi + u!2)T ~ 27rn, while Vn are somewhat analogous to 
"counter- propagating mode" resonances found in modu- 
lated traps [10|| and imply 2ujnT ~ 27r. Contrary to the 
suggestion of [8| where no Liapunov exponents were cal- 
culated, we find that none of the L{n, I) resonances have 
any |A| > 1. They are all stable, including L(l, 1), by far 
the strongest of all. But counter-intuitively, they are as- 
sociated with a much stronger BEC response, even after a 
very long-time (100 kicks) than the nonlinear resonances 
Nn and Vn which are unstable. 

The reason for this is clear from Figl2ljb). The mode 

2 populations from the CPE (for both A^i and N3) grow 
exponentially for a finite time, then decay exponentially; 
the inset shows this behavior on a log-scale. The map 
with BL corrections here is quantitative for only the first 
10-20 kicks, so we cannot model this behavior from iJ^^^ 
alone. But it is tempting to attribute it to regimes where 
either the A, A* or the 1/A, 1/A* eigenvalues are predom- 
inant. The cluster of interacting resonances N2,vi and 
V2 lies in a region of very low ground state depletion (it 
lies in the mimimum of the dominant L{1, 1) tail). The 
map ([U, corrected by BL terms H^^'' in Eq. ([4]) (with 
the lowest 7 modes), reproduces very well the charac- 
teristic Fano-like profiles seen in all three peaks, while 



the uncorrected map produces only symmetric resonance 
profiles. To our knowledge neither the exponential os- 
cillations, nor the Fano profiles are seen in comparable 
non-equilibrium BEC dynamics. While not fully under- 
stood, they indicate that the (5-kicked systems offer new 
and experimentally accessible BEC dynamics. 

While the Talbot time 17 = resonances have been 
proposed for metrological applications (e.g., for measure- 
ment of gravity), the similarly sharp BEC cut-off suggests 
analogous possibilities as it provides a sharp excitation 
threshold. One could envisage a set of ring-BECs subject 
to kicks; at threshold kicking frequencies, a very small 
perturbation would determine whether a particular ring 
is left in an ^ = state or acquires a large I = 1 com- 
ponent. Moreover, even under a very weak rotation, a 
BEC in a large ring would occupy a higher momentum 
(vortex) state I. Then the system could provide a sensi- 
tive probe of rotation for small changes in the resonance 
frequency between different angular momentum states. 



* Electronic address: 't.monteiro@ucl.ac.uk' 
[1] S. Fishman, I. Guarneri, and L. Rebuzzini, Phys. Rev. 

Lett 89, 0084101 (2002); S.Wimberger, LGuarneri and 

S.Fishman Phys. Rev. Lett 92, 084102 (2004). 
[2] M.K. Oberthaler et ai, Phys. Rev. Lett. 83, 4447 (1999); 

C. Ryu et ai, ihid. 96, 160403 (2006); M. Sadgrove et ai, 

ibtd. 99, 043002 (2007); G. Behinaein, et al, ibid. 97, 

244101 (2006); L Dana et ai, ibid. 100, 024103 (2008); 

J.F. Kanem et al, ibid. 98, 083004 (2007). 
[3] A. Buchleitner et al, Phys. Rev. Lett. 96, 164101 (2006); 

M. Saunders et al, Phys. Rev. A 76, 043415 (2007); M. 

Lepers V. Zehnl, and J-C Garreau, 77, 043628 (2008). 
[4] Z-Y Ma, K Burnett, MB d'Arcy, SA Gardiner Phys Rev 

A 73, 013401 (2006). 
[5] S. Wimberger, R. Mannella, O. Morsch, and E. Ari- 

mondo, Phys. Rev. Lett. 94, 130404 (2005). 
[6] C. Zhang, J. Liu, M.G. Raizen, and Q. Niu, Phys. Rev. 

Lett. 92, 054101 (2004). 
[7] L. Rebuzzini, R.Artuso, S.Fishman, LGuarneri, Phys. 

Rev. A 76, 031603(R) (2007). 
[8] J. Reslen, C. E. Creffield and T. S. Monteiro, Phys. Rev. 

A 77, 043621 (2008). 
[9] J.J. Garcia-Ripoll, V.M.Perez-Garcia and P Torres, 

Phys. Rev. Lett. 83, 1715 (1999); Yu Kagan and L A 

Manakova, cond-mat/0609159 (2006). 
[10] C.Tozzo, M.Kramer and F.Dalfovo, Phys.Rev.A 72 

023613 (2005); P.Engels, C.Atherton and M.A.Hoefer 

Phys. Rev. Lett. 98 095301 (2007). 
[11] D.M. Stamper-Kurn et al, Phys. Rev. Lett. 81, 500 

(1998). 

[12] E. Hodby et al, Phys. Rev. Lett. 86, 2196 (2001). 

[13] N.Katz, E.E.Rowen,R.Ozeri and N Davidson, Phys. Rev. 

Lett. 95, 220403 (2005) (2005). Phys. Rev. A 77, 033602 

(2008); E.E.Rowen, N.Bar-Gill and N.Davidson Phys. 

Rev. Lett. 101, 010404 (2008). 
[14] A. Sinatra, C. Lobo and Y. Castin J. Phys. B. 

At.Mol.Phys 35, 3599-3631 (2002). 



5 



[15] We introduce dimensionless units x = 2kLx',t = 8urt' , 
where the lattice recoil frequency oJr = hki/2M and 
g = (^^x- Since ftz,L = tt, L -» 27r and the BEC 
occupies < X < 27r. Also 2gn — > ^ and the Bogoli- 



ubov spectrum loI fc'^/2(fc^/2 + g/n). 
[16] K. Henderson et al, Europliys. Lett. 75, 392 (2006); 

Phys. Rev. Lett. 96, 150401 (2006). 
[17] A. Smerzi et al, Phys. Rev. Lett. 79, 4950 (1997). 



